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In models with several large extra dimensions and fundamental Planck scale of the order of 1 TeV, 
black holes can be produced in large numbers at LHC energies. We compute the charged hadron 
spectra obtained from the decay of black holes created in pp and Pb+Pb collisions at LHC. We 
show that hadrons from black hole decay dominate at transverse momenta pt ~ 30 — 100 GeV/c 
compared to usual QCD processes and black hole production signals are easy to identify in hadron 
transverse momentum spectra. Furthermore we show that a measurement of the charged hadron 
spectra probes Planck scales up to 5 TeV for any number of extra dimensions. 
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I. INTRODUCTION 



In scenarios with large extra dimensions the funda- 
mental Planck scale could be in the TeV range Jj. One 
of the most striking consequences of a low fundamen- 
tal Planck scale is the possibility of producing black 
holes and observing them in future colliders or in cos- 
mic rays/neutrino interactions I^I^Qil^. 

If gravity propagates in c? = 4 -t- n dimensions while the 
other fields are confined to a 3-brane, the 4-dimensional 



Planck scale is given by M|,; = Mp+^Vn = G:^l^Vn, 
where Mp is the fundamental Planck scale in 4 -I- n di- 
mensions and Vn = (27ri?)" is the volume of the n-torus 
that describes the compact space. For large size of the 
extra dimensions R, the fundamental scale Mp can be as 
low as in the TeV range. The existence of large compact 
dimensions leads to deviations from Newtonian gravity 
at distances of the order of R, as well as strong effects 
of Kaluza-Klein excitations of the graviton on various 
processes at high energies. These effects impose con- 
straints on the scale Mp, depending on the number of 
extra dimensions. For n < 4 the strongest limits are 
given by astrophysical processes like supernovae cooling 
and neutron star heating: Mp > 1500 TeV for n = 2 
and Mp > 100 TeV for n — 3. Cosmological considera- 
tions give similar bounds for n = 2 and n = 3 and they 
imply Mp > 1.5 TeV even for n = 4. These bounds con- 
tain, however, a larger degree of uncertainty and model- 
dependence. Non-observation of black hole production 
in cosmic ray interactions also imposes constraints for 
n > 4 at the level of Mp > 1 TeV. Present collider limits 
are typically below 1 TeV for any number of large extra 
dimensions. 

Given all these constraints, we will concentrate here 
on the scenario with n = 6 and Alp ~ 1 — 5 TeV. We 
will also discuss the dependence of our results on n and 
Mp. One thing to note is that the constraints mentioned 
above are derived for a toroidal compactification where 
all the large extra dimensions have the same radius. For 
different types of compactification the limits could actu- 
ally be considerably relaxed. 

Black hole production and evaporation can be de- 
scribed semiclassically and statistically when the mass 



of the black hole is very large compared to the funda- 
mental Planck mass. When the mass of the black hole 
approaches Mp one expects quantum gravity effects to 
become important. We want to explore only the param- 
eter space where the semi-classical treatment is justified. 
The total available energy at LHC is 14 TeV. Black holes 
with masses of this order can be produced in pp colli- 
sions. For Pb-I-Pb collisions, the black holes produced 
would have masses up to 5.5 TeV. These masses should 
be high enough above the Planck scales considered here 
for the semi-classical description to be valid Q . 

In Ref. , it was suggested that the geometrical cross- 
section would be exponentially suppressed. Detailed sub- 
sequent studies 't'I,^ did not confirm this proposal and 
showed that the geometrical cross-section is modified 
only by a numerical factor of order one. In Ref. , it was 
shown that even including the exponential suppression of 
the cross-section, the production rates are still high. We 
will use here the geometrical cross-section. 

These black holes decay very rapidly. The decay oc- 
curs in several stages. For the purpose of detecting 
black hole events, the most important phase is the semi- 
classical Hawking evaporation, since it provides a large 
multiplicity of particles and a characteristic black-body 
type spectrum. Most of this Hawking radiation is on 
the brane 0, producing all Standard Model particles. 
Because most of the Standard Model degrees of freedom 
come from strongly interacting particles (quarks and glu- 
ons), hadrons will be the dominant signal for the events 
where black holes are formed. 

In this paper we compute the transverse momentum 
distributions of charged hadrons at mid-rapidity obtained 
from the evaporation of black holes produced in pp col- 
lisions and Pb-|-Pb collisions at LHC energies. We show 
that in pp collisions the black hole events produce a 
large number of hadrons and dominate over the QCD 
background at transverse momenta above around 30-100 
GeV/c, where they can be clearly measured. The results 
have a weak dependence on the number of large extra di- 
mensions, but depend quite strongly on Mp. The signal 
is big enough to detect even for Mp ~ 5 TeV. For Pb+Pb 
collisions the energy available is lower, so one can only 
produce lower mass black holes and probe lower Plank 
scales. However, the rates could be higher due to the 
large number of binary collisions. We also discuss some 
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of the possible uncertainties that would affect our results. 

II. BLACK HOLE PRODUCTION AND DECAY 

In a high energy parton-parton collision, the impact 
parameter could be smaller than the Schwarzschild radius 
in d dimensions for a black hole with mass Mbh- 
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This leads to the formation of a semi-classical d- 
dimensional black hole of size rh much bigger than the 
fundamental Planck scale Mp, but much smaller than 
the size of the large extra dimensions R. 

With the above assumptions, the black hole production 
in a parton-parton collision is given by the geometrical 
cross-section cfbh — tj^ . Then the cross-section in a pp 
collision is obtained by folding in the parton densities: 
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where xi and X2 = M^h /{^is) are the momentum frac- 



tions of the initial partons and x 
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scale in the parton distribution functions f{x,Q ) is cho- 
sen to be = 1 /r^ . The results depend only weakly on 
the choice of this scale. We use CTEQ6M 10] for the 
parton distribution functions. 

The radiation rate into Standard Model particles is 
given by a thermal distribution in 4 dimensions: 
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with the black hole temperature: 
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where the sum is over all Standard Model particles and gi 
is a statistical factor, counting the number of degrees of 
freedom. The sign in the denominator is -I- for fermions 
and — for bosons, cr^^s are the gray body factors, which 
depend on the spin s of each particle. We first approx- 
imate these by ai^s — ^s^i, where Tg are constant [Tlj | 
(ri/2 = 2/3, Ti = 1/4, To = 1). A black hole acts as 
an absorber with a radius somewhat larger than r/j , such 
that Ak can be written as 0: 
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We compare these results with those obtained by using 
the recently computed gray body factors given in eq. (47) 
and (53) of Ref. ■ The differences in the final results 
are found to be around 30%. 

For the emission into gravitons, which are d- 
dimensional, the rate is given by: 



dE 
~dt 



1 



(27r)° 



Lugiffid k 



(7) 



where ai ^ A^. This will be much smaller than the 
emission rate into SM particles Q and we can neglect it 
when calculating the lifetime. 

The lifetime of the black hole is then obtained by in- 
tegrating eq. (PJ . Assuming no mass evolution during the 
decay, we get: 
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III. HADRON SPECTRUM 

We compute the cross section for inclusive charged 
hadron production from the partons produced in the 
black hole decay from 
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where Zmin = '^p/V^ smd z = p/pc and the decay distri- 
bution is: 



Ec 



dN _ 1 p'^u^jgcCTcTBH 

d^Pc ~ (27r)3 gP'aU^/TBH ± I ' 



(10) 



where 7 is the Lorentz gamma factor and u = 
(7,0,0,(2:1 — X2)y/s/{2M)) takes into account that the 
black hole is not produced at rest, but with a small ve- 
locity. Here p and E refer to hadrons, while Pc and E^ 
are for partons. 

We choose the scale of the fragmentation function 
Z?^(z,Qj) to be the final transverse momentum pT of 
the hadrons. The KKP fragmentation function 13] is 
used to get the final charged hadrons from the partons 
produced in the evaporation of the black hole. The KKP 
fragmentation function is only parametrized in the range 
of 0.1 < z < 1.0 and 1.4 < Q/ < 100 GeV. We need to 
access small values of the transverse momentum fraction 
z = Phadron/Pparton, as wch as large Qf, for the partons 
from black hole decay. For the large Qf, we explicitly 
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FIG. 1: Cross section for inclusive charged hadron production 
from 10- dimensional black holes in pp collision at -^s = 14 
TeV compared to the QCD background as a function of Mp. 
Black hole masses are integrated over 12TeV < Mbh < \/s. 
NLOpQCD predictions are plotted for the scales Q — Qf — 
2pt,Pt,Pt/2. 

evolve the KKP fragmentation function from the scale 
Qf — 100 GeV up to the desired values (in this case 
up to 10 TeV) using DGLAP equations For the 
z < 0.1 range, we use small-z fragmentation function by 
Fong and Webber [l5] which is based on the coherent 
parton branching formalism, which correctly takes into 
account the leading and ncxt-to-lcading soft gluon singu- 
larities, as well as the leading coUinear ones. It was found 
that the predicted energy dependence of the peak in the 
^ = ln(l/2;) distribution agrees very well with the e+e" 
annihilation data up to cm. energy of 200 GeV [l^ . 

Fig. n] shows the cross section for inclusive charged 
hadron production from black holes for several values of 
Mp ranging from 1 to 5 TeV in pp collision at LHC, com- 
pared to the expected spectrum of hadrons from QCD. 
LHC will be sensitive enough to detect the QCD hadrons 
up to pt around 400 GeV/c. The black hole signal is 
much bigger than the QCD one starting at px ^ 50 — 200 
GeV/c, depending on the Planck scale. It can be seen 
that even for Mp as high as 5 TeV there is a consider- 
able signal above background at pr ~ 200 GeV/c. At 
higher px the background is practically inexistent, while 
the black hole signal is still very large, as seen in Fig.|2] 

We show in Fig. [3 the cross-section for inclusive 
charged hadron production from black hole decay for 
Mp = 2 TeV, compared with the QCD background. It is 
reasonable to assume that black holes with masses only 
slightly higher than Alp can be produced. However, as 
previously discussed, our semiclassical description of the 
black holes is only valid for black hole masses much big- 
ger than the fundamental Plank scale. Integrating over 
the mass of the black hole starting at a low minimum 
value (close to Mp) would assume the validity of this 
treatment beyond its range of applicability. We choose 



FIG. 2: Cross section for inclusive charged hadron production 
from black holes in pp collision at -^5 = 14 TeV for n = 4 
and n — 6 and different minimum masses of the black holes 
produced. 

to integrate starting at a much higher black hole mass. 
Clearly, black holes with masses lower than our cut-off 
can be produced, but results obtained with a cut-off at 
Mp would have to take into account quantum gravity 
effects which are unknown. We show here results for 
M]^Jp = 10 TeV and M]^}}' = 12 TeV. It can be seen 
that including lower mass black holes gives considerably 
higher rates. Consequently, we consider our approach to 
be a 'conservative' one: our results are an underestimate 
of the actual signal and our qualitative conclusions al- 
ways hold, while the actual quantitative results could be 
much higher than our estimates, making the signal eas- 
ier to detect. Even for high mass black holes the signal 
clearly dominates over the QCD background in the re- 
gion above 100 GeV/c, where it can be easily seen in the 
experiments. Including lower mass black holes gives a 
bigger signal for all momenta and also drives the signal 
above the background even for lower momenta, of the 
order of tens of GeV. 

We also study the dependence of the results on the 
number of extra dimensions and show that it is very 
small, as can be seen in Fig|2 

The QCD background 17] is shown for different choices 
of the scale used in the structure and fragmentation func- 
tions (we use pt,Pt/2, 2pt)- The dependence on this 
scale is very weak in the high transverse momentum re- 
gion. For the black hole signal, the same change in Qf 
leads to differences of up to a factor of 2 in the results, 
which does not change any of our conclusions. 

We notice that even though there are significant 
changes in the overall rate of hadrons produced, the 
transverse momentum dependence of the hadrons does 
not change much when changing Mp or Mg^. This is 
not the case at the parton level. Changing Mp or il/]^^' 
the temperature of the black hole is modified and con- 
sequently the spectrum of the emitted particles is differ- 
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ent. Even though we can still see this for partons, the 
hadronization washes out most of the effect. We conclude 
that we cannot get a direct determination of the temper- 
ature of the black holes from the hadron spectrum. One 
could attempt to do that by looking at the spectrum 
of photons and electrons in the black hole event, which 
preserves the black body radiation type of spectrum 
but is considerably lower than the hadron signal because 
photons and electrons are only a small fraction of the 
particles produced in the black hole evaporation. In that 
case, one would be forced to consider black holes with 
lower masses in order to obtain a detectable signal. 

We conclude from here that black hole events will be 
easily detected in the hadron spectra at high px- The 
values of px for which the signal becomes higher than 
the background would give an indication on the values of 
Mp and Mbh that this signal corresponds to. 



IV. BLACK HOLES IN PB+PB COLLISIONS 

To compute the spectra in the case of Pb+Pb col- 
lisions, we multiply the expression in Eq. by the 
Glauber profile density Taa(6) = / d^rTAir)TA{\b - 
r|), where lA(r) = J dzp[r,z), normalized such that 
J (PrTjiir) — A, A being the atomic mass number and 
p(r, z) the nuclear density for which we take the Woods- 
Saxon distribution. This factor gives an enhancement 
for the production cross-section. For example, TAA{b = 
2fm) = 28 mb~^ in the case of Pb-|-Pb collisions at im- 
pact parameter b — 2 fm. However, y/sJjN = 5.5 TeV 
in this case, so only lower mass black holes can be pro- 
duced and a smaller parameter space can be probed. We 
do not include shadowing effects, since black holes dom- 
inate at high X, where these are negligible. In addition, 
it is demonstrated that nuclear modification of parton 
distribution is getting smaller when we go to larger scale 
from Q2 ^ 2.25 GeV^ up to = (lOOGeV)^ We 
have confirmed that this also holds true for much higher 
scales up to Q = 30 TeV by evolving the nuclear parton 
density in 'l&l and conclude they are negligible. 

In Fig. Owe show the results for Pb-|-Pb collisions for 
Mp = l TeV and M]^jp = 5 TeV. For these parameters, 
the signal is still significant and can be easily detected. 
Due to the fact that is only 5.5 TeV, the higher 

scales are no longer accessible in this type of experiment. 

In Pb-|-Pb collisions, the black hole is expected to 
be produced in a dense medium of quarks and gluons, 
therefore we need to take into account the interactions 
of the partons produced in the decay of the black hole 
with the quarks and gluons around it, for example, as in 
Ref. [3 . The energy loss is expected to be small at high 
transverse momentum. At LHC energies, for px ^ 100 
GeV/c, where the black hole signal clearly dominates 
over the background, the energy loss was found to be 
small, (about 5% effect) 

However, energy loss has significant effects at pt below 
10 GeV in the QCD spectrum at LHC For the 




Pj (GeV/c) 

FIG. 3: Transverse momentum distribution for charged parti- 
cle at mid-rapidity from black hole decay in Pb-(-Pb collision 
at impact parameter 6 = 2 fm at LHC. NLOpQCD calculation 
with the scale Q = Qf = pr is also shown. 



hadrons coming from black holes we also expect the effect 
to be small in the high momentum region. However, there 
is a possibility to have enhanced particle yield around 
Pt ^ 10 GeV/c, because the hadron spectra is much 
fiattcr than that of the QCD spectra and feedback from 
the emitted gluons could be non-negligible. If this is 
the case, the black hole signal could be also identified in 
the lower pT region, in addition to the high px one. It 
could happen that the signal becomes higher than the 
background at these values of pr- Even if the signal is 
somewhat lower than the background, it is still large, 
such that the experiment would detect a large number of 
hadrons from signal+background, even at pt's of tens of 
GeV. 

We do not include interactions of the black hole itself 
with the surrounding particles. One can imagine taking 
into account possible absorption of these particles by the 
black hole, which would affect the decay of the black 
hole. This is an interesting issue and is presently under 
investigation 



V. CONCLUSIONS 

In summary, we have computed the transverse momen- 
tum spectra for high px charged hadrons from decay of 
black holes produced in pp collisions at y/s — 14 TeV as 
well as in central Pb+Pb collisions at y^sJ^N = 5.5 TeV. 
We have shown that the hadrons from black holes are 
detectable and dominate the background for px above 
about 100 GeV for fundamental Planck scales up to 5 
TeV. Our results are conservative, as they only take into 
account very high mass black holes. Including black holes 
with lower masses gives even stronger signals. The value 
of Pt at which the signal becomes bigger than the back- 
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ground is determined by Mp and Mbh considered. 

We have neglected the evolution of the mass of the 
black hole during the decay. If we take into account 
that, as the mass decreases, the temperature increases, 
we would get a somewhat harder spectrum. In the same 
time, the lifetime would decrease compared to our esti- 
mate, so that our curves would move slightly down and 
to the right. However, all the qualitative features previ- 
ously discussed will remain the same. The last stage of 
the decay, when the mass of the black hole has decreased 
to almost Mp is not understood, since it requires a full 
quantum gravity description. This is why a full consid- 
eration of the mass evolution during decay is not really 
possible. We have not taken into account the angular 
momentum of the black hole and the phase of the decay 
when the black hole would just loose this angular mo- 
mentum. Also, there are a few additional particles pro- 
duced in the initial stage of the black hole decay, when 
the black hole looses the quantum numbers of the par- 
tons that produced it. In ;8j it has been shown that 
classical gravitational radiation could be important and 
consequently the actual black hole mass is smaller than 
that of the center of mass energy of the parton collision. 
This would reduce the number of very high mass black 
holes produced. However, even in case of a small fraction 
of initial energy going into black hole production, as long 
as Mbh ^ Mp, there is still a large, observable signal 
since production of black holes with small mass is large. 

If there are additional degrees of freedom around 100 
GeV (new particles), so below TpH, they should be taken 



into account and they would lead to a small decrease in 
the lifetime. Their contribution would be slightly sup- 
pressed, just as for the top quark, Ws, Z's and Higgs, 
for which the masses are no longer much smaller than the 
temperature. These issues introduce some uncertainty in 
the numerical results, but our main conclusions are un- 
affected. 

We would like to note that the QCD background we 
show in the graphs is computed at y = 0. For high ra- 
pidity this background is actually much smaller, while 
the black hole signal is the same for all rapidities. This 
would indicate that by looking in the high rapidity region 
one would enhance the signal to background ratio even 
further. 

In conclusion, we have shown that the charged hadron 
distribution in pp and Pb-|-Pb collisions at LHC energies 
provides an unique probe of black hole production and 
the physics of extra dimensions for Planck scale up to 5 
TeV and for any number of extra dimensions. 



VI. ACKNOWLEDGMENTS 

We thank Radu Roiban for useful discussions. We 
would also like to thank P. Aurenche and J. P. Guil- 
let for providing us the code for calculating NLOpQCD 
hadron distributions. This work has been supported in 
part by the DOE under Contracts DE-FG02-95ER40906 
and DE-FG03-93ER40792. 



[1] N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. 
Lett. B 429, 263 (1998); I. Antoniadis, N. Arkani-Hamed, 
S. Dimopoulos and G. R. Dvali, Phys. Lett. B 436, 257 
(1998). 

[2] S. B. Giddings and S. Thomas, Phys. Rev. D65, 056010 
(2002); S. Dimopoulos and G. Landsberg,Phys. Rev. 
Lett. 87, 161602 (2001). 

[3] T. G. Rizzo, JHEP 0202, Oil (2002). 

[4] P. C. Argyres, S. Dimopoulos, J. M. March-Russell, Phys. 
Rev. Lett. B441, 96 (1998); T. Banks, W. Fischler, 
hep-th/9906038 S. Nussinov, R Shrock, Phys. Rev. D59, 
105002 (1999). 

[5] J. L. Feng and A. D. Shapere, Phys. Rev. Lett. 88, 
021303 (2002); L. Anchordoqui and H. Goldberg, Phys. 
Rev. D 65, 047502 (2002); R. Emparan, M. Masip and 
R. Rattazzi, Phys. Rev. D 65, 064023 (2002); S. Iyer 
Dutta, M. H. Reno, I. Sarcevic, Phys. Rev. D 66, 
033002 (2002); L. Anchordoqui, J L. Feng, H. Gold- 
berg, A. D. Shapere, Phys. Rev. D 65, 124027 (2002); 
A. Chamblin, G. C. Nayak, Phys. Rev. D66, 091901 
(2002), L. Anchordoqui, H. Goldberg, hep-ph/0209337 
K. Cheung, Phys. Rev. Lett. 88, 221602 (2002), Phys. 
Rev. D66, 036007 (2002); T. Han, G. D. Kribs, B. McEl- 
rath, hep-ph/0207003 M. Cavaglia, hep-ph/0210296 
A. Ringwald and H. Tu, Phys. Lett. B 525, 135 (2002); 
M. Kowalski, A. Ringwald and H. Tu, Phys. Lett. B 529, 



1 (2002); S. Hossenfelder, S. Hofmann, M. Bleicher, 
H. Stocker, Phys. Rev. D66, 101502 (2002); R. Casa- 
dio, B. Harms, Phys. Lett. B487, 209 (2000), Phys. 
Rev. D64, 024016 (2001), Int. J. Mod. Phys. A17, 4635 
(2002). 

[6] M. B. Voloshin, Phys. Lett. B518, 137 (2001); 
M. B. Voloshin, Phys. Lett. B524, 376 (2001). 

[7] D. M. Eardley, S. B. Giddings, Phys. Rev. D66, 044011 
(2002); H. Yoshino, Y. Nambu, Phys. Rev. D66, 065004 
(2002); S. D. H. Hsu, hep-ph/0203154 S. N. Solo- 
dukkhin, Phs. Lett. B533, 153 (2002), S. Bilke, E. Li- 
partia, M. Maul, hep-ph/0204040 A. Jewicki, J. Thaler, 
Phys. Rev. D66, 024041 (2002). 

[8] H. Yoshino, Y. Nambu, Phys. Rev. D67, 024009 (2003). 

[9] R. Emparan, G. T. Horowitz, R. C. Myers, Phy. Rev. 
Lett. 85, 499 (2000). 
[10] J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, P. Nadol- 
sky, and W. K. Tung [CTEQ Collaboration], JHEP 0207 
(2002) 012. 

[11] D. N. Page, Phys.Rev.D 1 3, 198 (1976) 

[12] P. Kanti, J. March-Russell, |hep-ph/0212199| Phys. Rev. 

D66, 024023 (2002). 
[13] B. A. Kniehl, G. Kramer and B. Potter, Nucl. Phys. B 

582, 514 (2000); Nucl. Phys. B 597, 337 (2001). 
[14] We used the code Qcdnum version 16.12 for the 

DGLAP evolutions which can be obtained from 



6 



http://www.nikhef.nl~h24/qcdnum/ 

[15] C. P. Fong and B. R. Webber, Nucl. Phys. B355, 54 
(1991), V. A. Khoze and W. Ochs, Int. J. Mod. Pys. 
A12, 2949 (1997). 

[16] W. Braunschweig et al. [TASSO Collaboration], Z. Phys. 
C 47, 187 (1990); D. Akrawy et al. [OPAL collaboration], 
Phys. Lett. B247, 617 (1990); P. Abreu et al. [DELPHI 
Collaboration], Eur. Phys. J. C 18, 203 (2000) [Erratum- 
ibid. C 25, 493 (2002)] 

[17] P. Aurenche, M. Fontannaz, J. Ph. Guillet, B. A. Kniehl, 
E. Pilon, M. Werlen, Eur. Phys. J C9, 107 (1999); P. Au- 
renche, M. Fontannaz, J. Ph. Guillet, B. A. Kniehl, 
M. Werlen, Eur. Phys. J C 13, 347 (2000); 



[18] K. J. Eskola, V. J. Kolhinen, P. V. Ruuskanen, Nucl. 
Phys. B535, 351 (1998); K. J. Eskola, V. J. Kolhinen, 
C. ,A. Salgado, Eur. Phys. J. C. 9, 61 (1999). 

[19] M. Gyulassy, I. Vitev, and X. N. Wang, Phys. Rev. Lett. 
86, 2537 (2001); S. Jeon, J. Jahhan-Marian, and I. Sarce- 
vic, nucl-th/0208012 hep-ph/0207120 T. Hirano and 
Y. Nara, Phys. Rev. C 66, 041901(R) (2002); T. Hirano 
and Y. Nara, nucl-th/0211096 

[20] 1. Vitev and M. Gyulassy, Phys. Rev. Lett. 89, 252301 
(2002) 1. Vitev, hep-ph/0212109 

[21] 1. Mocioiu, Y. Nara, I. Sarcevic, in preparation. 



